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Topological insulators are noninteracting, gapped fermionic systems which have gapless boundary 
excitations. They are characterized by topological invariants, which can be written in many different 
ways, including in terms of Green's functions. Here we show that the existence of the edge states 
directly follows from the existence of the topological invariant written in terms of the Green's 
functions, for all ten classes of topological insulators in all spatial dimensions. We also show that 
the resulting edge states are characterized by their own topological invariant, whose value is equal 
to the topological invariant of the bulk insulator. This can be used to test whether a given model 
Hamiltonian can describe an edge of a topological insulator. Finally, we observe that the results 
discussed here apply equally well to interacting topological insulators, with certain modifications. 
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I. INTRODUCTION 

The time-reversal invariant topological insulators of 
recent interest, and the integer quantum Hall systems 
of longstanding interest (for recent reviews at a variety 
of levels, please see Ref. [T]), are now known to be just 
two elements in a classification table of all noninteract- 
ing fermionic systems,'^ which identifies ten symmetry 
classes of topological insulators. All of these systems are 
gapped in the bulk and possess a bulk topological invari- 
ant. They are also supposed to have topologically pro- 
tected gapless excitations at the boundary. These gapless 
excitations are often taken as the most significant, even 
the defining, property of topological insulators. 

A variety of arguments can be given supporting the 
existence of gapless edge states based on bulk properties, 
some more detailed (such as Laughlin's original argument 
designed for the integer quantum Hall effect!^ and some 
more qualitative. For example, consider the case of a 
model of noninteracting fermions that supports both a 
topologically trivial phase and a nontrivial phase, with 
the phase transition driven by varying a parameter in the 
system's Hamiltonian. At the transition, corresponding 
to a special value of the parameter, the excitation gap 
must close. A boundary between two samples in the two 
phases can then be seen as a domain wall across which 
the parameter varies spatially through its special value. 
This "spatial phase transition" should then also result in 
gapless excitations, which form in the vicinity of the spe- 
cial value. Such arguments, while correct, tell us nothing 
about what kind of gapless edge excitations might form 
at the boundary of the topological insulator. 

Here we rely on the method of Green's functions to 
give a general, quantitative argument that proves the ex- 
istence of the edge states for all topological insulators. 
Moreover, our argument shows that the edge states are 
described by their own topological invariant, analogous 
to the winding number of a vortex in a superfluid or to 
the charge of a particle as measured by Gauss' law (and 
similar to the defect invariants of Teo and Kan^, whose 
value must be equal to the invariant of the bulk insula- 



tor. This edge invariant vanishes if the edge states are 
gapped, giving a simple way to see why they must be 
gapless. The existence of this invariant gives us a tool 
to test whether a particular Hamiltonian can describe an 
edge theory, as it is straightforwardly calculable for any 
Hamiltonian under consideration. 

Before we derive this argument later in this article, 
let us describe the edge topological invariant. A d- 
dimensional topological insulator that is translationally 
invariant, so that the d-dimensional momentum p'' is a 
good quantum number, possesses a d-dimensional (bulk) 
topological invariant iV^.I^If the topological insulator has 
an edge, translation invariance in the direction perpen- 
dicular to the edge is lost, and the good quantum number 
is the (d— l)-dimensional momentum p''"^ parallel to the 
edge. Let us take one of the components of this momen- 
tum, say Prf-i, and fix it at some large value A. Although 
we expect the edge to have gapless excitations at some 
momentum (and so is not an insulator), the Hamiltonian 
at fixed pd-i = A is gapped if A is sufficiently large (as 
a function of the remaining d — 2 momenta), and so de- 
scribes a (d— 2)-dimensional insulator. We will show that 
this insulator is a topological insulator, with the topolog- 
ical invariant Nd-2{^)- Finally, we will show that 

iVrf = iVrf_2(A)-7Vd_2(-A), (1) 

which constitutes the main result of this article. 

As an example, consider the three-dimensional edge of 
a four-dimensional time-reversal invariant topological in- 
sulator with spin-orbit coupling. Suppose the excitations 
localized at the boundary are described by the Hamilto- 
nian 

'H = v ^ PaO-Q - li, (2) 

a—x.y,z 

where CTq are the Pauli matrices and v is Fermi velocity, 
which satisfies the appropriate symmetry (time-reversal 
invariance) H{p) = UyH* {—p)ay. Its single-particle ex- 
citation spectrum is 

e±{^) ^±vJpI+pI+pI- (3) 
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and its zero energy excitations occur on a sphere of radius 
IJl/v, as shown in Fig. [I] Let us see that this is indeed 
the edge of a topological insulator. We fix = ±A to 
obtain the two-dimensional Hamiltonian 

H± = V {px(7x + Py(^y ± Acr^) - /i, (4) 

again as illustrated in Fig. [T] This Hamiltonian, un- 




FIG. 1. The Fermi surface of radius /x/w of the topological 
insulator discussed in the text, and the surfaces Pz = ±A on 
which the edge topological invariant is computed. 

derstood as describing a two-dimensional system, is well 
known. For example, it has been studied as a model of 
quantum Hall transitions.'^ The two-dimensional topo- 
logical invariant of this Hamiltonian is known to be equal 
to its Hall conductance in units of /h (note that Eq. Q, 
unlike Eq. ([2]), breaks time reversal invariance), and, as- 
suming A > /i/f , evaluates to 



A^2(A)-7V2(-A)-1. 



(5) 



Therefore we conclude that this H \s a. valid model of 
the boundary of a four-dimensional topological insulator. 
This illustrates the usefulness of Eq. ([T ) . 

This example assumes an unphysically large bulk di- 
mension, and the reader may doubt that it has any rele- 
vance for real- world physics. However, it is knowrP that 
the experimentally realized time-reversal-invariant topo- 
logical insulators in 2 and 3 dimensions can be viewed 
as descending from a four-dimensional parent, simply by 
setting one or two of the momenta to zero. Similarly, the 
three-dimensional boundary just discussed can be seen as 
the parent of the boundary theories of those real systems. 
In particular, it follows from Eq. ([2| that the boundary 
theory for the two dimensional time-reversal invariant 
topological insulator is H — v cr^Px — fJ-^ while that of 
the three dimensional time-reversal invariant insulator is 
H = V {(7xPx + '^yPy) ~ We will return to this "dimen- 
sional reduction" later in this article. 

The Green's function formalism we use is quite pow- 
erful: it remains applicable even when interactions are 
added to the system, since Green's functions exist for in- 
teracting as well as for noninteracting systems. Here one 



must comment that the power of topological invariants 
often lies in the fact that they represent the response of 
noninteracting systems to external electromagnetic field. 
Once the interactions are turned on, the topological in- 
variant written in terms of interacting Green's functions 
may or may not represent the response, and this could be 
considered an impediment to using the Green's function 
formalism to describe interacting topological insulators. 
However, the main result, Eq. ([!]), is correct regardless 
of the presence of interactions as will be clear from its 
derivation. If one takes as the most physically relevant 
property of topological insulators is the existence of edge 
states, then since the relation Eq. ([I]) relates the topo- 
logical invariant to the property of the edge, this rela- 
tion, and not the response, can in principle be taken as a 
starting point of the application of this formalism to the 
interacting topological insulators. 

Yet the interpretation of Eq. ( ij may change once the 
interactions are turned on. As discussed in a recent pa- 
per by one of us,'^ in the presence of interactions the 
(eigenvalues of the) Green's functions may have zeroes 
localized at the edge, in addition to the poles that indi- 
cate the usual edge states. From the point of view of the 
topological invariants, zeroes are similar to poles, and can 
result in a nonzero edge invariant even in the absence of 
zero-energy edge states.fl^- In fact, it can be shown that 
the recent result of Ref. [TTl where a model topological 
insulator loses its edge states in the presence of inter- 
actions, is due to the replacement of the edge states by 
the zeros of the Green's function. The study of the ef- 
fects of interactions using the formalism developed here 
appears to be a useful direction of further research, but 
goes beyond the scope of this article. In what follows 
we mostly assume that the fermions do not interact, al- 
though throughout this article we point out what exactly 
changes if interactions are taken into account. 

To connect this article to previous publications, we re- 
mark that the relation Eq. ([T]) was discussed for d — 2 
in Ref. [T31 Subsequently, it was extended for d = 1 and 
discussed in the presence of interactions in Ref. M Here 
we extend it for all topological insulators, in any number 
of dimensions. 

We proceed as follows: in section |ll] we derive Eq. ^ 
for systems with no symmetry other than translational 
(class A in the Altand-Zirnbauer classificatiorP^HH) ^ 
which only have integer-valued topological invariants in 
even spacial dimensions. Then in section [HT] we describe 
the effects of discrete symmetries on the invariants; and 
then describe how these results imply the relation be- 
tween bulk and surface physics for Z2 topological sys- 
tems. Finally in section [TV] we discuss systems with chi- 
ral symmetry. Appendix |A| contains the derivation of a 
key step in our derivation of Eq. ([T]), and Appendix [B| 
presents a more complete discussion of ID topological 
insulators. 
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II. TOPOLOGICAL INSULATORS WITHOUT 
ANY SYMMETRIES 

We start with a translationaUy invariant topological in- 
sulator with a single particle Green's function Gafiii^, p)- 
If interactions are absent, the insulator can be described 
by the Hamiltonian 



(6) 



Here dj^p creates a fermion with momentum p in d dimen- 
sions, and the species label a runs over the spin, bands, 
as well as particle/hole space for a superconductor (sum- 
mation on indices should be always understood) . In that 
case, the Green's function is simply given by 



(7) 



Once interactions are switched on, simple expressions of 
this kind are no longer available, but Green's function 
can still be defined in a standard way (see, for example. 
Ml- 

The bulk topological invariant is 



j diod^ptlG-^da^G . . .G-^da,G, (8) 



where a runs over (w, p) , and e is fully antisymmetric 
tensor, with ecjpi...pi — +1 (the other nonvanishing com- 
ponents are obtained by permutation). The constant Cd 
is given by"^ 



C, = -(27rz)-W2)-i(d/2)!/(d+l)! 



(9) 



The spatial dimension d is even, as Nd vanishes when d 
is odd (by antisymmetry of e and cyclicity of the trace 
tr); we will later introduce a different expression, valid 
for chiral systems in odd dimensions, that we denote by 
the same symbol. 

This topological invariant always evaluates to an inte- 
ger, and gives the Hall conductance in two-dimensional 
spacep^ at least in the absence of interactions — it re- 
mains a topological quantity (an integer) even with in- 
teractions. This is true simply because it measures the 
winding of the map from a (d -f l)-dimensional space 
(cj, p'') to a space of matrices G; it is known that such a 
map can be topologically nontrivial if d is even. More- 
over, it is clear that in the absence of interactions, to 
change the value topological invariant one needs to de- 
form the Green's function in such a way that either G or 
G^^ becomes singular. For noninteracting systems, this 
is only possible if w = and % has zero energy states as 
follows from Eq. ([7|; this results in infinite G at cj = 
and the closing of the gap in the spectrum. For inter- 
acting systems G itself can acquire zero eigenvalues at 
Lu — 0,- SO that G^^ is singular; thus interacting systems 
can change the value of their invariant Eq. ([s]) without 
ever closing the gap. All of this matches known results. 



Now we introduce some technology, inspired by 
Volovik.EsI Consider a d-dimensional system with a do- 
main wall of dimension d — 1 ; the Hamiltonian varies in 
the direction perpendicular to the domain wall in such a 
way that far from the domain wall the Hamiltonian de- 
scribes a topological insulator with the invariant or 
iVj", on the right or left side of the domain wall (what 
this means precisely will be defined below). Since there 
are only d— 1 good momenta, it is described by the mixed 
Green's function G(a;, p''"^; s, s'), where s is the coordi- 
nate normal to the boundary and is effectively another 
matrix index. A Fourier transform with respect to s — s' 
produces the Wigner (-Weyl) transform G(a;, p*^, s)j ^^ l ^^ l 
with s = (s + s')/2 and pd the momentum conjugate 
to s — s' . Far from the domain wall, this will become 
s-independent and coincide with the bulk Green's func- 
tion (that is, with translation invariance the Wigner and 
Fourier transforms are the same). 

Importantly, two objects exist which can be interpreted 
as the inverse of G. The "true" inverse satisfies 



ds'K^fiis, s') Gp^is', s") = 5^^6{s - s") (10) 



(in this expression the lo and p"* ^ dependence of K and 
G was suppressed for brevity). For noninteracting sys- 
tems 



K{uj, p''"^ s, s') =iuj- H{iu, p'*"^ s, s') 



(11) 



but for interacting systems one should rely solely on 
Eq. ( 10 ) to calculate K. At the same time, matrix (or 
local Jmverse G~^ satisfies 

G-i(u;,p'',s)G^^(c^,p^s)=<5„^. (12) 

With translation invariance along s, K — G^^, but in the 
presence of a domain wall K ^ G~^ [here K{uj, p'', s) is 
the Wigner transform of K{uj, p'*"^; s, s')]. 

With these tools, we define the (d + 2)-dimensional 
vector 

Gd tiG-^da,G...G-'da,G (13) 

in the space (w,p'',s). Remarkably, the divergence of 
this vector is zero except where G (or G~^ if interactions 
are present) is singular, as can be checked by direct dif- 
ferentiation. This is completely analogous to the electric 
field in electrostatics with point charges or charged sur- 
faces. Here, the sources of n are singularities (or zeros) 
of G. The analogy to electrostatics is productive; we 
will essentially measure the charge of a singularity with 
Gauss' law, integrating the flux of n over a surface that 
surrounds the charge. 

For the rest of this section, for simplicity we will always 
refer to these as "singularities" having in mind that in the 
absence of interactions G cannot have zeros. We will also 
keep in mind that if interactions are present, singularities 
will imply either infinite or zero G. 
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The singularities of G may occur at multiple points of 
the d + 2 dimensional space or on surfaces in the space 
spanned by (w, p'', s). We denote these fi. 

Because G'(a;,p'',s) reduces to the bulk Green's func- 
tion far from the boundary (domain wall) , the definitions 
Eqs. Q and ( 13 1 mean that we can compute Nd as the 
flux of n through the surface s = i for large L. That 
is, the difference between topological invariants on either 
side (right and left ^ say) of the boundary, 



= j dujdf^p[ns{s ^ L) 



n,{s^~L% (14) 



is just the flux of n through the combined surface 
(a;,p'^,±L), shown in Fig. 2. On the other hand, and 




FIG. 2. The surfaces a — ±L on which A^^ is computed as 
the flux of n. 



this is the crucial point, because dana = we can use 
any surface we like to compute — so long as it 
encloses the singular surfaces of G. Denoting by Sf. any 
(d + 1) -dimensional "Gaussian" surface (for example, a 
sphere) surrounding fi we flnd 



N 



(15) 



Eq. (15) constitutes the flrst half of our argument.^ 



do] 



To relate the flux of n to the edge properties of the 
system, we construct a d-dimensional vector r out of the 
mixed Green's function G(w, p'^"^; s, s'), which allows us 
to define a topological invariant that directly captures 
the behavior of the edge states. We define 



— Cd-2 Cao. 



,Tr/^oa,„Go. 



■oKoda,_,G. (16) 



The convolution {AoB){s,s") = Jds' A{s, s')B{s' , s") 
that appears here is really a generalization of matrix mul- 
tiplication to functions of two coordinates s, s' . That is, 
the coordinates s, s' are treated as matrix indices of the 
functions A and B, and G; of course, these functions also 
carry the ordinary matrix indices a, /3 that label fcrmion 
species, etc, and which are being summed over in the 
usual way. Then Tt A = J ds tr A{s, s) is just the trace 



for these generalized matrices. The index a here runs 
over (a;,p''~^), the space of the vector r. 

The vector r, just like the vector n above, is divergence- 
free, except at singularities of the mixed Green's function 
G. These singularities can also be interpreted as "elec- 
trostatic" sources for the vector r emanating from them. 

It is clear that the singular surfaces fi — the sources 
of n — become the sources of r as well when projected 
from the {d + 2)-dimensional space (w,p'',s) onto the 
d-dimensional edge space (a;,p'^~^). Since the Green's 
function G — [iuj — H] ^ can only be singular where w — 
and H has zero energy eigenstates, these sources are 
confined to w = and form surfaces (or points) in the (d— 
l)-dimensional space spanned by p'*"^, which we identify 
with the edge Fermi surfaces and Dirac points and denote 
Fi (having in mind that in the presence of interactions 
surfaces and points of zero G play an equivalent role). 

Remarkably, we can prove the following statement: 



(17) 



Here Sp^ is the (d — 2)-dimensional surface surrounding 
the Fermi surface (or Dirac point) Fi. The proof of this 
statement is given in Appendix [A] and is quite involved 
(although for d = 2 the proof is significantly simpler; it 
is given in Ref. [TS]) . It relies on the approximation of 
the domain wall being smooth. However, the corrections 
to this equation form a series expansion in powers of the 
gradient dgG. Since both sides of this equation are inte- 



gers due to Eq. (151, small corrections to it must vanish. 



It follows from Eqs. (15 1 and (17) that 



(18) 



In words, this says that the Fermi surfaces (and Dirac 
points) on the edge are characterized by "topological 
charges," or fluxes of r emanating from these surfaces 
and points. This charge is equal to the difference in the 
bulk topological invariants on either side of the bound- 
ary. Eq. ( 18 ) is the quantitative statement of the fact 



that gapless excitations are required at a boundary be- 
tween bulk insulators with different values of the topo- 
logical invariant Nd (in the presence of interactions, there 
can be zeroes insteacP). Indeed, in the absence of such 
excitations, r is divergence-free everywhere and its flux 
is always zero. 

What remains is to interpret this charge as a difference 
of topological invariants, as discussed in the beginning 
of this article. So long as the Fermi surface does not 
traverse the Brillouin zone in the Pd-i direction, we can 
choose Sp to be the surfaces (w, p''"^, ±A) for suitable 
A, instead of the spheres closely surrounding the Fermi 
surface(s) that we assumed until this point [see Fig. |3]. 

The flux through Sp is then the difference of the fluxes 
of r through these two surfaces. In turn, those fluxes can 
be reinterpreted as the difference Nd-2{-^) — A^d-2(— A), 
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d-2 




Pd-1 



FIG. 3. The surfaces p"" 
as the fiux of r. 



±A on which Nd-2 is computed 



where 



Nd-2{Pd-l 



did d 



d-2 



(19) 



is a (d — 2)-dimensional topological invariant calculated 
in the space (w, p''"^) with pd-i fixed, as can be verified 
with Eq. ( jlG] ). The main result of this article Eq. ([T]) 
(where is assumed to be for simplicity) immediately 
follows. 

In the simplest case d — 2 this result was derived in 
Ref. HSJ where iVo(A) — A'o(— A) was interpreted as a 
number of energy levels crossing zero as the momentum 
along the one-dimensional edge changes from —A to A, 
corresponding to the standard picture of edge states in 
the integer quantum Hall effect. In higher number of 
even dimensions d, this result's interpretation is more 
complicated, but the result is nonetheless useful, as in 
the example given in the Introduction. 

This concludes the discussion of topological insulators 
without symmetries (systems of class A). 



III. TOPOLOGICAL INSULATORS WITH 
TIME-REVERSAL OR PARTICLE-HOLE 
SYMMETRY 

Now let us apply this formalism to systems in the nine 
remaining classes of topological insulators with symme- 
try, beginning from non-chiral insulators. Nonchiral insu- 
lators, termed classes AI, All, C and D, are those which 
possess either time-reversal or particle-hole symmetry (it 
can be helpful to refer to the "periodic table" of topolog- 
ical insulators in Ref s. [Hand [3] for the discussion of this 
section and the next). We have seen that with no symme- 
try (class A), bulk insulators are classified by an integer 
in even dimensions. The same invariants describe the 
other nonchiral classes in even dimensions. However the 
presence of discrete symmetries can force them to vanish 



in some dimensions, as follows. Classes AI and All repre- 
sent time-reversal (T) invariant systems with integer and 
half-integer spin respectively, and the Green's function 
satisfies the constraint 



(20) 



where Ut is a unitary matrix such that U^Ut = ct = +1 
for AI and = — 1 for All {G'^ is the transposed Green's 
function). The same constraint holds for G. Putting this 
into Nd, relabeling p — p, taking the transpose, using 
cyclicity, and relabeling indices, we see that this forces 
Nd = when d ^ An + 2 (an anticyclic permutation is 
odd in d = 4n -I- 2). Similarly, classes C and D have 
particle-hole (C) symmetry 



G(c^, p, rs) = -Ul G^i-Lo, -p, ,s) Uc 



(21) 



where UqUc = ec = +1 for class D and ec ~ — ^ for class 
C. An additional minus sign in the integrals is generated 
by w — > —uj, so that is nonzero only in spaces of 
dimension d — An + 2. All of this matches what was 
established elsewhere using a different language's! 

The conclusion here is that if we are in a spatial di- 
mension where the topological invariant can be nonzero, 
the bulk-boundary correspondence Eq. ([T]) applies just 
as it does for systems without any symmetry. One might 
worry that the relation Eq. ([T]) involves not only but 
also Nd~2'i as we just saw, in systems with symmetries 
if Nd is nonzero, Nd-2 appears to be zero. However, in 
Eq. ([ij Nd-2 is calculated with p"^^^ fixed at some value 



A, so the Green's function no longer satisfies Eq. (20 1 or 
( [2T| ). Thus both sides of Eq. ([T]) can be nonzero, as they 
should. 

Let us now account for the appearance of Z2 invariants, 
which in time-reversal-invariant insulators which belong 
to class AH appear for example in d = 2 and d — 3. 
To explain these, we appeal to the dimensional reduction 
picture of Qi, Hughes, and Zhang,!^ in which the physi- 
cal system is considered as embedded within a space of 
higher dimension, with one or two of the momenta being 
fictitious additional parameters in the Green's function. 
The real physical system corresponds to the fictitious mo- 
menta set to zero. Under T (or C) almost every point in 
the momentum space has an image at the opposite mo- 
mentum, the exceptions being the time-reversal-invariant 
(TRI) points. As a result, if a Green's function has a 
point singularity at a non-TRI point, this contributes 2 
to Eq. (18), since each point and its image contribute. 



So, if Nd — N^ is odd and all singularities are points, 
some of the singular points must be located at the TRI 
points in the momentum space. Finally, it is possible to 
make sure, by choosing the appropriate extension of G to 
the unphysical momenta, that all singular TRI points are 
in the physical space. This leads to the conclusion that if 
the topological invariant in the extended space is an odd 
integer, there must be singular points, or as we saw edge 
states, in the physical (reduced-dimension) theory. 

The same arguments apply if, instead of a Fermi point, 
there is a Fermi surface surrounding a TRI point. The 
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case of a surface is also crucial for showing that gener- 
ally the dimensionality cannot be reduced by more than 
2. Indeed, generalizing a result due to HofavapS generic 
topological insulators characterized by integer invariants 
have two-dimensional Fermi surfaces at the boundary. A 
Fermi surface (spherical, for simplicity) centered on zero 
is defined by the equation pi +pl+pt = Pf^ with the 
remaining momenta, if any, being arbitrary. This can be 
dimensionally reduced by one or two by setting pz = 
or Py = Pz — 0. However, to reduce the dimensionality 
by more than 2 requires setting Px — Py — Pz — 0, which 
eliminates gapless excitations. 

An example of two-dimensional Fermi surfaces located 
at the three-dimensional edge of a four-dimensional topo- 
logical insulator with time-reversal invariance is shown 
in Fig. [4j Three Fermi surfaces are depicted, giving 



\ 



Py 



FIG. 4. Example of an edge Brillouin zone containing a 
time-reversal-invariant, multicomponent Fermi surface. The 
two smaller spheres each contribute to the three-dimensional 
edge topological invariant, but they do not intersect the time- 
reversal invariant planes of the Brillouin zone and so do not 
contribute to dimensionally reduced edges. 

Ni = iV2(A) - 7V2(-A) = 3. At the same time, two of 
those Fermi surfaces located at opposite momenta play 
no role if dimensions are reduced by setting, for example, 
Pz = 0. The remaining Fermi surface centered on the ori- 
gin still yields zero energy excitations even if = 0, thus 
in this example the edge state survives in lower dimen- 
sions. On the other hand, if 7V4 were even, then the cen- 
tral Fermi surface would necessarily be absent (or there 
would be a pair that could be deformed off the pz — 
plane), and there would be no edge excitations in lower 
dimensions. Here, as everywhere else, in the presence of 
interactions, the Fermi surfaces could be the surfaces of 
zero energy excitations as well as surfaces of zeroes of 
Green's functions. 

Note that if one breaks time reversal invariance by 
adding appropriate terms in the Hamiltonian, the exci- 
tations no longer have to be symmetric under reflection 
of momenta. Then the central Fermi surface can move 
off the center of the Brillouin zone and no longer con- 
tribute to the lower dimensional edge excitations. This 



is the mechanism by which breaking time reversal invari- 
ance removes edge states in lower dimensions (but not 
in the original 3-dimensional edge where edge excitations 
survive even if time reversal invariance is broken) . 

Finally, the dimensional reduction gives a trivial result 
in some dimensions because the bulk invariant Nd takes 
only even values for the dimensionally extended system, 
which washes out the Z2 structure. 

These arguments conclude our discussions of topolog- 
ical insulators in classes AI, All, C and D. 



IV. TOPOLOGICAL INSULATORS WITH 
CHIRAL SYMMETRY 

We also want to show that analogous results hold for 
the classes with chiral, or sublattice {S) , symmetry ( AIII, 
BDI, CII, CI and Dili). In these systems there is a matrix 
S such that 

G(w,p,s) = -EG(-a;,p,s)E, = 1. (22) 

The bulk invariant for chiral systems, analogous to 
Eq. ([8|, can be written as 



-ea,...aJd''ptTj:G-^da,G . . .G-^da,G, (23) 



where G is evaluated at cj = 0, a runs over the compo- 
nents of p, and the bulk dimension d is now odd.[23l Unlike 
the even-dimensional topological invariant Eq. (|8|, this 
expression does not treat momenta and frequency in a 
symmetric fashion and in fact is derived directly from the 
topological invariant written in terms of the Hamiltonian 
in Ref. Elby replacing V. — > G^^\^^q- An expression for 
Nd with d odd which involves integration over frequency 
and momenta is also possible,!^ however for the purpose 
of this article it is not needed. 

Expressions strictly analogous to Eqs. (13), (15), and 
( 16 ) can be formulated for the chiral case, leading to 



a relation just like the bulk-boundary correspondence 
Eq. (18). In particular. 



G-'da,G, (24) 



where a runs over (p'^,s), is a, {d + l)-dimensional 
divergence-free vector, while 

Tao ^^^ea,...a,_,TTj:Koda,Go---oKda,_,G (25) 

(a runs over p''^^) is a (d— l)-dimensional divergence-free 
vector. Here everything is evaluated at w = 0. 

Class AIII, without T or C symmetry, is potentially 
nontrivial in any odd dimension. Classes BDI, CII, CI, 
and Dili are characterized by the presence of both T 
and G symmetries, with er^c — V equal to -1-1 in the 
first two classes and —1 in the other two. The invariant 
Eq. (23) vanishes if ?/ = 1 and d = 3 + An, or if 77 = 
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— 1 and d = 1 + An, because E and G anticommute at 
a; = as a consequence of S'-symmetry and Ut'SUI, = 
rjYp" . Arguments similar to the ones employed above in 
nonchiral classes can now be used to see the emergence 
of the chiral Z2 invariant in appropriate dimensions. 

For example if d = 1, the vector r reduces to a scalar 
r = — TrE. As is well known, zero-energy states are 
eigenstates of E, with the eigenvalue +1 (right zero 
modes) or —1 (left). TrE, the difference of the number 
of right and left eigenstates, is then equal to the differ- 
ence of the bulk invariants on both sides of this zero- 
dimensional boundary. A reader may worry that the for- 
malism reported earlier in this article becomes somewhat 
degenerate at c? = 1, therefore we present a slightly more 
detailed derivation of this result in Appendix [B| 

This concludes the derivation of Eq. ([T]) for all classes 
of topological insulators in all spacial dimensions. 



V. CONCLUSIONS 

In this article we presented a derivation of the rela- 
tionship Eq. ([1]) between the bulk and the boundary of 
topological insulators. Even when there are no interac- 
tions and no disorder, this relation is quite useful and 
allows to test whether a particular proposed boundary 
theory can indeed be at the edge of a topological in- 
sulator. When interactions are present, the relationship 
Eq. remains true and still allows to relate the the bulk 
to the edge, although the edge states may get replaced by 
zeroes.^ In addition to the extension to interacting sys- 
tems already described, the method described here can 
also be extended to disordered topological insulators, by 
periodically repeating the finite size system. The results 
obtained in this way will be reported elsewhere. 



transform 



A{p^,s)= j dre^P^^A[s+-rs--). (Al) 

The technical tool we need in the following calculation is 
the Moyal product expansion, which allows us to write 
the Wigner transform of C = A o B (this stands for 
C{s, s") = / ds' A{s, s')B{s', s")) as 



C^AB + ^^e^,d^Ad,B 



(A2) 



where /i, stand for s, p^. In particular, this means that 
1 



K = G ' + -e,,.XpX.G i+. 



X^ = G-'d^G. (A3) 
This can be checked by demanding that K o G = 1, ap- 



plying Wigner transform and using Eq. (A2). 



The Moyal product expansion is a formal expansion in 
powers of the gradient of G, and the omitted terms in- 
volve higher gradients of G. As mentioned in the text, 
this does not limit the applicability of our results to slow 
domain walls; since we are deriving a relation between 
integer-valued quantities, any corrections from higher 
powers of the gradient must vanish (unless they also eval- 
uate to integers, in which case they would indicate extra 
topological structure) . 

We would like to relate the vector r to the vector 
n. The former involves Green's functions in real space, 
while the latter involves Wigner-transformed Green's 
functions. Therefore, we need to apply the gradient ex- 
pansion to 



Gd-2^~ 



First, consider the trace of a product: 



(A4) 
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Appendix A: Derivation of Eq. (17| 



TiAoB 



dsds' tr A{s,s')B{s\s) 



dsds' dpddp'^ 



(27r)2 
dsdpd 



d g~iPdis~s')-ip'^(s'-s) 



A I s + s'\ f , s + s' 
tr A [ Pd, ] B [pa, — — 



2tt 



tr A{pd,s)B{pd,s). 



(A5) 



This states that the trace is invariant upon changing to 
the Wigner- Weyl basis. The expansion we need is there- 
fore 



TrKoda.G^ 



'Koda,_,G 



dsdpd 
2tt 



tr KM, (A6) 



Given a matrix A{s, s') which is also a function of two 
variables s and s', one can introduce its Wigner (or Weyl) 



where M is the Wigner transform of 9^1 G°K ■ ■ ■od^^_^Go 

Koda,_A 
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M = 



2i 



2i 



2j ^tid-2,ad-2^^'^d-2 



K 



2i 



f^d-2 '^d-2 



dad-iG 



2i 



2i 



dad-^GK 



2i ' 



-'Hd-2,ad-2^'-'i^d-2\^^ y ~r 2j 



Eq. ( A7 1 comes from applying the Moyal product expan- 



sion, Eq. (A2), to every o-product. 



The notation is quite unwieldy, so we make the follow- 
ing notational simplifications. First, we suppress all e^j,; 
since we will only be keeping terms of first order in the 
gradient, there should be no confusion about how indices 



are contracted. Second, the factor e,. 



_j is suppressed 



as well, and every other symbol that carries an index is 
written in bold face, with the index suppressed: Xai — > X 
and — > d. It must be remembered that these symbols 



daa-iG - 



(A7) 



I 

form a completely antisymmetric tensor. In more for- 
mal language, these are one-forms. We can also rewrite 
products (and o products) 



(A8) 



Finally, we will write d—1 — n (which is odd by assump- 
tion) . 

Then the quantity to be computed, up to the overall 
constant, is 



TrX"= /^trX 
2tx 



AGK + j^d^dGd,{K •) + ^dGd^Kd, 



dG 



The bracket expands to 

{dGKY'-^ + - ^(dGG-i)"-2^^ [a^dGa^(G-i •) -f dG^^G^^a^] (dGG'^)^ 



3=0 



where K G ^ in the sum because we keep only first order in the gradient. Then 
" dsdpd 



TrX" = 



2n 

n-2 



tr j(ifdG)" 

- ^(G-MG)"-2-^G-i [df,dGd^{G-^ •) + dG^^G^^a^] dG(G-MG)^ 
1 " 



dsdpd 
2tt 

n-2 



1 " ' 

+ ^ E X""'"' [G-'d^dGd,{x'+') - XX^.G-^^,{Gx'+')\ 



3=0 



(A9) 



(AlO) 



(All) 



where we have expanded K in the first term according produces what is often termed a "weak" topological in- 
variant, that is, a topological invariant of lower dimension 
than the bulk system; we will ignore this term from now 



to Eq. (|A3l) and replaced x ^ G^MG and d^G'^ 

-Xti-G- 



Note that the the zeroth-order term in the expansion, on. 

Because n is odd, no sign is picked up under cyclic 
relabeling of the indices hidden in the bold-face notation. 



dsdpd 
2-K 



tr X" 



(A12) 



so we can use cyclicity of the trace. Therefore, we can 
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combine terms to obtain the first order quantity 



The combination in square brackets is simply dx^ so, 
shifting the summation index, 



n n — 2 

T ^ ^ tr x''-^x,XuX' + E X""'"' 

n-2 



(dx^)a.(x^). (A13) 



The factor dx^ contains double derivatives of G, which 
do not appear in our desired result. Therefore, we will 
separate out total derivatives in d, which produce no 
flux on the (closed) Gaussian surfaces Sp on which the 
topological invariant is built [see Eq. (17)]. To do this, 
note that d^ = and that, when using the product rule 
to expand these derivatives, a sign appears every time d 
moves past a %. This gives 

j'=0 

which implies 







J even 



jodd 



, (A14) 



tr x'^^'-'i'ixMxn 



since n is odd. Then 



dtrx" ^ '^Xt^d^ix^) jeven 
^-dtr x"-'-'x^dAxn - tr x""-' X^d^ix') + x'^-'-^^.d^X'^^) J odd 

I 

gives 



(A15) 



T = tr 



x'\,Xu + dE(-l)'tr x''-^-\^d,{x') 



T = Di + tr x"XmX. 

n-l j 
j — 1 m— 1 

= Di +tr x^'x^Xu 



E x''''-'x,d,{x'+') - tr x'^-'x^d^ix')] 



ji-i 

E 



j odd 



E(-i)^' 



tr x"-™XpX™-^9.X 



= Di+trx"XMX.+ E tr x"-"XmX'""'9^'X- 



m^2 
m even 



dE(-l)'tr x"-^-1xm5.(x^) +tr x"XmX. 
i=i 

+ E(-I)'trx"-^XM5.(X^'), (A16) With = dx. + XX. - X.X, 



Di+trx"XMX.+ E tr x"~"XMX'""'dx. 



m=2 
m even 



where we have shifted the summation index in the last 
term. The total derivative is of no further interest, so 
we set Di = dX:;=i (-Iptr x"-^-'Xp9,(x^). In the 
last term, expanding the derivative di, and using cyclicity 



E [trX^-^X^X^X.-trx^-^+'x^X^-'x.] ■ 

(A17) 



m^2 
m even 



10 



Now, because n is odd and m is even, 

m=2 
m even 

n-1 



m^2 
m even 



-trx"-"(dxp)x'"-'x. 

+ tr x"-"XmX"x. + tr x"-"XMX'"-'dx.] 

n-l 



m^2 
m even 



tr x"-"+'XmX™-'x.] 



(A18) 



where the second equahty follows from antisymmetry in 
fi and v after a relabeling of the summation index in the 
second term of the first equality. This means 

T = Di + tr x"XmX. 



^ 71—1 



tr x"-'"XmX""'x. 



m^2 
m even 



tr ^"-"XmX^X. 

71— m+1 m — 1 



- tr X'^^^+'X^X^-'X.}. (A19) 

We can collect the total derivatives as D2, shift the sum- 
mation index in the last term, and rewrite 2 tr x"X/jXi/ — 
tr x"XmXi^ - tr XmX"Xi^ to obtain 

1 " 

T = D2 + - ^ (-1)™ tr ^"-"XpX^X.- (A20) 



Including D2 with the suppressed terms and restoring 
e^^, this means that 



TrX" 



dsdpd 



X^(-ir6^,tr x"-"XmX™X. 



m=0 



(A21) 

This expression is sufficiently simple that we can return 
to the original, more explict notation; 



dsdpd 



^liu^a^ai ...an 



X tr Xai ■ • ■ XQ„-™XMXa„_„+i ■ • • Xa„Xi^ + 



(A22) 



The right-hand side has 2n\{n -|- 1) = 2(n -I- 1)! terms 
and is totally antisymmetric, while using the Levi-Civita 
symbol gives {n + 2)1 terms. Therefore, 



Tr O ■ ■■ O Xa 



1 r 

T^-i-i) J ^^^^'^ ^o-oai...a„+2^^ Xai ■ ■ ■ Xa„+2 



(A23) 



assuming that we choose to order the new indices appro- 
priately. 

Finally, it is important that the constant be correct: 



[d + iy. 



{2iii 



d/2 



d{d+l) 



-{2TTi) 



-# W2-1)! 
[d-l)\ 



47ri((i+ 1) 
so (substituting back d = n + 1) 



(A24) 



Cd-2eaoai...ad_iTrXai O • • • O Xa^_i 
= -;^—^J^^^Cd-2jdsdpdeaoai...aa+,tTXm ■ ■ ■ Xa^+i 

+ ... 

= Cd J dsdpd eooai...orf+i tr Xai • • • Xaa+i + (A25) 



or 



dsdpd Ua 



(A26) 



This is valid for values of a in (a;,p''~^), but does not 
make sense for a = s or a = pd since r does not have 
such components. 

We would like to calculate 



+1 



/ 



n. 



(A27) 



We deform the sphere Sf to the space 5*^ x M where Sp 
is the sphere in the space (a;,p''~^) and M is the entire 
space (M^) spanned by pd and s. This can be termed 
a "hypercylinder," since a cylinder is the product of a 
circle (analogous to our sphere Sp) and a straight line 
(analogous to the flat space M = M^). 

The flux of n through such a surface can be computed 
with only the components of n for which the relation 



Eq. (A26) holds. Part of that flux involves integration 
over dsdpd, which then naturally leads to 

JdSj+^-n^ JdSp^- Jdsdpdn = J dSj-^-r, (A28) 



where Eq. (A26) was used. This is the result asserted in 



the text as Eq. (17) 



Appendix B: Chiral systems in one dimensional 
space 



We start by writing Tr S as 

Tr S = Tr E 7^ o G. 



(Bl) 
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Notice that this is true at arbitrary w. We then take 



advantage of Eqs. (A3) and (A5) to rewrite this as 



— trE / dsdp(G-^d-sGG-^d„G ~ G-^dMG'^dsG) . 
Airi J ^ f f / 

(B2) 

The expression to be integrated is a total derivative and 
results, upon integrating, in 



TrE = iVi(s = -L) - iVi(s = X), 



(B3) 



where A^i, defined in Eq. (23), is evaluated at s = —L 
and s = L. This is indeed what is claimed in section ITVl 
What remains is to remark that while these expressions 
are written at arbitrary uj, the limit w — >■ is conve- 
nient as only in this limit Tr S counts zero energy states 
(and zeroes of Green's functions if there are interactions) 
localized at the boundary. This derivation can be used 
instead of the much more involved procedure reported in 
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